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A dynamical systems approach to competition of Saffman- Taylor fingers in a Hele-Shaw channel 
is developed. This is based on the global study of the phase space structure of the low- dimensional 
ODE's defined by the classes of exact solutions of the problem without surface tension. Some 
simple examples are studied in detail. A general proof of the existence of finite-time singularities 
for broad classes of solutions is given. Solutions leading to finite-time interface pinch-off are also 
identified. The existence of a continuum of multifinger fixed points and its dynamical implications 
are discussed. The main conclusion is that exact zero-surface tension solutions taken in a global 
sense as families of trajectories in phase space spanning a sufficiently large set of initial conditions, 
are unphysical because the multifinger fixed points are nonhyperbolic, and an unfolding of them 
does not exist within the same class of solutions. Hyperbolicity (saddle-point structure) of the 
multifinger fixed points is argued to be essential to the physically correct qualitative description of 
finger competition. The restoring of hyperbolicity by surface tension is discussed as the key point 
for a generic Dynamical Solvability Scenario which is proposed for a general context of interfacial 
pattern selection. 



I. INTRODUCTION 

The Saffman-Taylor (ST) problem Jl] § has played a 
central role for several decades as a prototype system in 
the study of interfacial pattern formation |q-|il L partic- 
ularly concerning the issue of pattern selection [ l|,^2|dl4[ . 
Despite its elongated existence, the problem continues to 
pose new challenges with the focus now on its dynam- 
ical aspects fLq] . In this sense, the ST problem is be- 
coming instrumental once more in gaining insights into 
possibly generic behavior, due to its relative simplicity 
in the context of morphologically unstable interfaces in 
nonequilibrium systems. 

Full understanding of the analytical mechanisms lead- 
ing to steady state selection by surface tension as a sin- 
gular perturbation in the problem was not completely 
achieved until the late eighties Jl6|-p3| and the result- 
ing scenario, usually referred to as Microscopic Solv- 
ability (MS) |i|-||, has currently become a paradigm for 
many other systems for instance in free dendritic growth 
P, p!0p^| ]. Such solvability analysis, however, is strictly 
static, in the sense that it is concerned with the existence 
and linear stability of stationary solutions. The impor- 
tance of dynamics in the process of selection was pointed 
out in Refs. [f25H27| where it was argued that the Saffman- 
Taylor finger solution was not the universal attractor of 
the problem if the displacing fluid has a non-negligible 
viscosity. More recently, the traditional MS scenario of 
selection has not been exempt of some controversy in 
connection with the dynamics of the zero surface tension 
problem |28| p4| . The singular effects of surface tension 
on the dynamics have been pointed out as a rather subtle 
and challenging issue [pa p7l and the possibility of some 



extension of the MS scenario of selection to the dynamics 
has been suggested || j^p8|]3fj| l . In any case, the study 
of the dynamics of morphologically unstable interfaces in 
the context of Laplacian growth or, more generally, of 
diffusion-limited growth of interfaces in nonequilibrium 
conditions, has been rather elusive to analytical treat- 
ment due to the highly nonlinear and nonlocal character 
of the equations. Simplified models of Laplacian growth 
in the absence of surface tension have thus been studied 
in some detail Efl-p2], and very recently a new impulse 
to the problem of Laplacian growth without surface ten- 
sion, with focus on DLA-like self-similar growth, is de- 
veloping after new analytical insight s |43| - [45| , leading to 
somewhat controversial conclusions [|I6[. For the viscous 
fingering problem, however, the basic question remains 
as to what extent the zero surface tension problem does 
capture the physics of the fingering dynamics. The direct 
motivation of this question is the existence of large classes 
of explicit time dependent solutions of the zero surface 
tension which are nonsingular and span a great variety of 
possible interface morphologies. To what extent these so- 
lutions are qualitatively or quantitatively describing the 
dynamical behavior or real systems is an outstanding, 
nontrivial question. As we will see, in some cases they 
may evolve to the correct asymptotic solution with the 
wrong dynamics. In others, despite the fact that the solu- 
tions are unstable to arbitrary perturbations, they may 
accurately describe the time evolution of systems with 
nonzero surface tension. The difficulty of the question 
lies to a large extent in the precise way to formulate it. 

The present paper expands and elaborates in depth 
several aspects which were first pointed out in Ref . |38| . 
Our central objective is to contribute to elucidate the role 



of surface tension in the dynamics of finger competition 
of the ST problem in the channel geometry, beyond the 
obvious role of keeping the interface stable and free of sin- 
gularities and the selection of a stationary state. That is 
to elucidate the qualitative differences between the non- 
singular solutions of the zero surface tension problem and 
the real problem with a possibly small but nonzero sur- 
face tension. Our main contribution is the development 
of a new approach to the issue in terms of the ideas and 
concepts of the theory of Dynamical Systems (DS). With 
this general point of view, we will study in detail some 
specific classes of solutions of the zero surface tension 
problem. As a byproduct, this will provide some inter- 
esting results in the context of Laplacian growth concern- 
ing the interplay of finger widths and screening effects. 
However the most important issue will be the discussion 
of the physical deficiencies of the nonsingular solutions of 
the zero surface tension problem. This will be addressed 
with the aim at the maximum possible generality. As we 
will see, the comparison of the problem with and without 
surface tension is necessarily qualitative in nature, so it 
is important to pose questions in a framework which is 
at the same time qualitative and mathematically precise. 
Such framework is the theory of Dynamical Systems. The 
use of this conceptual tool will help us formulate precise 
questions to which we can give an answer. From the 
above results and within this spirit we will reformulate 
the issue of a possible extension to dynamics of the MS 
scenario of steady state selection, and suggest a possible 
answer to that. 

Common understanding of the finger competition pro- 
cess (sometimes referred to as finger coalescence) leading 
to the selected steady state is usually based on qualita- 
tive screening arguments |38| . In some cases these have 
been shown to be too naive, in accord to the recent find- 
ings of stationary finger solutions with coexisting unequal 
fingers |39|] , as we will discuss in more detail below. In 
a first effort to develop precise conceptual tools for a 
quantitative and qualitative characterization of finger- 
ing dynamics, a topological approach was developed in 
Refs. [ [25|]27| to address the nontrivial effects of viscosity 
contrast in the dynamics of finger competition and their 
long time asymptotics. The basic insight was to focus on 
topological properties of the physical velocity field in the 
bulk, quantified by the existence of topological defects. 
In this paper we will develop a different global viewpoint, 
which focuses on topological properties of the flow in the 
(infinite-dimensional) phase space of interface configura- 
tions, rather than the (two-dimensional) velocity field. 

In this paper we will analyze in detail some exact so- 
lutions of the Saffman- Taylor problem with zero surface 
tension. Hereinafter we will refer to this case as the ide- 
alized problem, as opposed to the regularized one, which 
will denote the problem with a small but finite surface 
tension. Is is known that the idealized ST problem is ill- 
posed as an initial- value problem |35|]. Nevertheless, one 
of the crucial facts that makes the ST problem attractive 
from an analytical point of view is the existence of rather 



broad classes of exact time-dependent solutions of the 
problem without surface tension J47|-[52| . Some classes of 
those solutions are known to develop finite-time singu- 
larities in the form of cusps and are thus of no interest 
to the physics of viscous fingering, but still a remark- 
ably broad class of solutions is free of singularities and 
therefore physically acceptable, in principle. The basic 
question is then what would be the effect of a small but fi- 
nite surface tension to those solutions. This question was 
first raised in Ref. |33| where it was shown that for some 
classes of initial conditions, the effect of surface tension as 
a perturbation could be considered as basically regular, 
while for other initial conditions the singular character of 
the perturbation showed up dramatically in the dynam- 
ics. In other configurations, such as for circular geometry, 
surface tension has also been shown to behave as a regular 
perturbation [p4[] . Indeed, in view of the morphological 
diversity which is included in the known nonsingular so- 
lutions, one may be tempted to believe that, since such 
solutions remain smooth for all the time evolution, they 
should stay close to the solutions of the regularized prob- 
lem as do — * 0. Siegel and Tanveer |36|,p7f have shown 
that, contrarily to what happens in other more familiar 
singular perturbations in physics, in the case of Hele- 
Shaw flows that is not the case, and, in general, the ide- 
alized and the regularized solutions differ from each other 
at order one time. In the remarkable contribution of Rcfs. 
|36| , p7| , however, only simple examples of single-finger 
evolutions are considered, so the extent to which those 
conclusions can be extrapolated to multifinger configura- 
tions still requires a careful analysis pq). Furthermore, 
even though the idealized and the regularized solutions 
differ significantly after a time of order unity which is ba- 
sically independent of surface tension, one could still ar- 
gue that the qualitative evolution is basically unaffected 
by surface tension if the finger width is not too differ- 
ent from the selected one in the regularized case. There- 
fore, the possibility that some classes of solutions or some 
particular dynamic mechanisms are basically insensitive 
to surface tension remains open. This question directly 
motivated the work of Ref. p8| which first made use of 
a dynamical systems approach and which we will here 
pursue further. The crucial aspect to be exploited is the 
fact that the integrable classes of initial conditions define 
finite-dimensional invariant manifolds of the full (infinite- 
dimensional) problem, so it makes sense to study the re- 
sulting low-dimensional dynamical systems and compare 
them with properly defined finite-dimensional subsets of 
the regularized problem. With this analysis we will clar- 
ify in what precise sense the nonsingular exact solutions 
of the idealized ST problem are, in general, unphysical. 
Once settled the unphysical nature of a broad class of 
solutions, a natural question to address in whether a se- 
lection principle is associated to the surface tension regu- 
larization, which can be understood as a dynamical gen- 
eralization of the MS scenario. We will address this point 
in the light of our results and discuss how and in what 
sense such dynamical MS could be formulated. 



The rest of the paper is organized as follows. In Sect. 
U the equations describing Hele-Shaw flows in channel 
geometry are recalled, together with the conformal map- 
ping formulation used to study the problem. The fin- 
ger competition phenomenon is described and the rel- 
evant quantities to characterize finger competition are 
presented. In Sect. Ill the dynamical systems approach 
to the problem is introduced. In Sect. |V| the minimal 
class presented in Ref. p8[ is revisited. The dynamics of 
this class, namely, its phase portrait is studied in detail. 
The comparison of this phase portrait with the known 
topology of the physical problem reveals that the mini- 
mal class dynamics is unphysical. The main reason for 
this unphysical behavior is the existence of a continuum 
of fixed points. In Sect. [V| and [VJ various generaliza- 
tions of the minimal class are introduced. We pay a spe- 
cial attention to the perturbation of the minimal class 
that removes the continuum of fixed points but keeps 
the dimensionality of the phase space unchanged. These 
solutions contain the main general features of zero sur- 
face tension, and it is shown that they do not describe 
in general the correct dynamics for finite surface tension. 
The main reason is that the dynamical system that they 
define lacks the saddle-point structure of the multifinger 
fixed point, necessary to account for the observed finite 
surface tension dynamics. In Sect. VII we discuss the 
precise role of zero surface tension solutions and their rel- 
evance to an understanding of the dynamics of Hele-Shaw 
flows. In this section a Dynamical Solvability Scenario is 
proposed and discussed as a generalization of MS theory. 
Finally, in Sect. VIII we summarize our main results and 
conclusions. 



II. FORMULATION OF THE PROBLEM AND 

CHARACTERIZATION OF FINGER 

COMPETITION 
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where T means that the quantity is evaluated on the in- 
terface, v n is the normal component of the velocity of the 
interface, k is the curvature, n is the unit vector normal 
to the interface and do is a dimensionlcss surface tension 
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defined by d = 12 a y W i ; where T4o is the fluid veloc- 
ity at infinity. Eq. (0) is the Laplace pressure jump at 
the interface due to local equilibrium, written in terms 
of the velocity potential. Here we consider the case with 
the non- viscous fluid displacing the viscous one, therefore 
the pressure in the non-viscous fluid is constant, and we 
choose it to be zero. Eq. (0) is the continuity condition, 
that states that the interface follows the motion of the 
fluid. We assume periodicity at the sidewalks of the chan- 
nel. Except for configurations symmetric with respect to 
the center axis of the channel, periodic boundary con- 
ditions define different dynamics from the more physi- 
cal case of rigid sidewalk (with no- flux through them). 
Strictly speaking our case describes an infinite periodic 
array of unit channels. We will argue that nothing es- 
sential is lost with respect to competition in a rigid-wall 
channel, while the analysis is significantly simplified. 

We use conformal mapping techniques to formulate the 
problem [Q. We define a function f(w,t) that confor- 
mally maps the interior of the unit circle in the complex 
plane u> into the viscous fluid in the physical plane z = 
x + iy. We assume an infinite channel in the x direction. 
The mapping f(u>,t) must satisfy d u f(oJ,t) ^= inside the 
unit circle, |w| < 1, and moreover, h(u>,t) — f(uo,t)+hxu) 
must be analytic in the interior. We define the complex 
potential as the analytic function $ = ip + itp, where the 
harmonic conjugate ij> of <p is the stream function. The 
width of the channel is W — 2ir and the velocity of the 
fluid at infinity is Voo = 1. It can be shown that the 
evolution equation for the mapping f(ui,t) reads 



Consider a Hele-Shaw cell of width W in the y- 
direction and infinite length in the x-direction, with a 
small gap b between the plates. The fluid flow in this 
system is effectively two dimensional and the velocity v 
obeys Darcy's law 
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where p is the fluid pressure and /i is the viscosity. We 
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define a velocity potential ip — ~-^jjP, and assuming that 
the fluid is incompressible (V • v = 0) we obtain the bulk 
equation to be the Laplace equation 



vV = o. 



(2) 



This must be supplemented with the two boundary con- 
ditions 



tp\r = d K\ r 



(3) 



dtf(w,t)=wd u f(fa,t)A 
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where A[g] is an integral operator that acts on a real 
function g(4>) defined on the unit circle. A[g] has the 
form 
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and on the unit circle u> 



it reads 



A[s]L =e i*=s($ + iH^] 



(6) 
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where H<p[g] is the so-called Hilbert transform of g(<ft) 
defined by 



*M = ^- p I \(0)cot(^-)d9 (8) 
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where P stands for the principal value prescription. The 
complex potential $ satisfies 

$(w,t) = -lnu; + doA[K] (9) 

where the curvature k given in terms of f(e 1 ^, t) is 
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The evolution equation (pj) written on the unit circle 



lj = e l * is 
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where f((f>, t) = f(e z ^, t). In the zero surface tension case 
do = the integro-differential equation (O) reduces to a 
much simpler equation, and the evolution of f(4>,t) for 
do = is then given by 



Re{ify/(<M)&/*(& *)} = !■ 



(12) 



The direct motivation of the present study is that, de- 
spite the fact that neglecting surface tension is in princi- 
ple incorrect from a physical standpoint, the do = case 
can be solved explicitly in many cases p] JIq , pO| including 
solutions which, although being unstable, they exhibit a 
smooth and physically acceptable (nonsingular) behav- 
ior, quite similar to what is observed in experiments and 
simulations of the full problem. 

Before proceeding to the description of the general ap- 
proach and its application to specific solutions, let us first 
introduce some ideas and definitions which will be help- 
ful in further discussion. To quantify finger competition 
it is useful to define individual growth rates of fingers 
p8[ . In simple situations like those considered in the pa- 
per, the growth rate of a finger can be simply defined (in 
the reference frame moving with the mean interface po- 
sition) as the peak-to-peak difference of the stream func- 
tion between the maximum and the minimum which are 
adjacent to the zero of the stream function located at (or 
near) the finger tip ||27[ (the definition can be generalized 
to more complicated situations). According to this def- 
inition, one assigns a nonzero growth rate to the finger 
if the tip advances at a velocity which is larger than the 
mean interface position. Looking at individual growth 
rates one can easily distinguish two different stages of 
the dynamics in the process of finger competition. A first 
stage characterized by the monotonous growth of all in- 
dividual finger growth rates and a second one dominated 
by the redistribution of the total growth rate among the 
fingers. We call these two stages growth and competition 
regimes respectively. For a configuration of two differ- 
ent fingers, which is practically the only one addressed 
throughout this paper, during the growth regime the two 
fingers develop from small bumps of the initially flat in- 
terface, while the total growth rate Ai/V(£), defined as 
AipT(t) = Atpi(t) + Aifj2(t), grows until it reaches a value 
close to its asymptotic one AVr(oo). The decrease of the 



growth rate of one of the fingers signals the outcome of 
the competition regime: there is a redistribution of flux 
from one finger to the other one. We also define the 
existence of successful competition as the ability to com- 
pletely suppress the growth rate of one finger. A finger is 
dynamically suppressed of the competition process when 
its growth rate Aip is reduced to zero. 



III. DYNAMICAL SYSTEMS APPROACH 



The theory of Dynamical Systems is a mathematical 
discipline which studies ordinary differential equations 
or flows (and also difference equations or maps) with 
stress on geometrical and topological properties of so- 
lutions [p6[. The approach is sometimes referred to as 
qualitative theory of differential equations. The focus 
is not on the study of individual solutions or trajecto- 
ries of the differential equation, but on global properties 
of families of solutions. This point of view has become 
very fruitful in searching for universality in the context 
of nonlinear phenomena. 

A dynamical systems approach seems thus appropriate 
to study in a mathematically precise way, the qualita- 
tive properties of the dynamics of our problem, and the 
qualitative differences associated to the presence or ab- 
sence of surface tension. One of the important concepts 
in dynamical systems theory is that of structural stabil- 
ity, which captures the physically reasonable requirement 
of robustness of the mathematical description to slight 
changes in the equations. Roughly speaking, a system 
is said to be structurally stable if slight perturbations of 
the equations yield a topologically equivalent phase space 
flow. Although the structural stability 'dogma' must be 
taken with some caution |36| , a structurally unstable de- 
scription of a physical problem must be seen in principle 
as suspect. When a dynamical system (DS) depends on 
a set of parameters, the bifurcation set is defined as those 
points in parameter space where it is structurally unsta- 
ble. In this case the structural instability at an isolated 
point in parameter space is the property necessary for the 
system to change its qualitative behavior. At a bifurca- 
tion point, adding perturbations to the equations to make 
the system structurally stable is called an unfolding H . 
For dimension higher than two, the mathematical defini- 
tion of structural stability is usually too stringent. For 
the purposes of the present discussion and most physical 
applications it is sufficient to consider the notion of hyper- 
bolicity of fixed points, which in 2 dimensions is directly 
associated to structural stability through the Peixoto the- 
orem |56| . A fixed point is hyperbolic when the linearized 
flow has no marginal directions, that is, all eigenvalues of 
the linearized dynamics are nonzero. We will see that the 
non-hyperbolicity of the double- finger fixed point (in gen- 
eral the n-equal-finger fixed point) and the non-existence 
of an unfolding of it within the known class of solutions 
is at the heart of the unphysical nature of this class of 



solutions. 

In the approach to the Saffman- Taylor problem with 
concepts of DS theory, there is, however, an impor- 
tant additional difficulty in the fact that our problem 
is infinite-dimensional and unbounded. In similar spa- 
tially extended systems, such as described by PDE's, it 
is customary to project the dynamics onto effective low- 
dimensional dynamical systems based on the so-called 
center manifold reduction theorem [ p6[ . This is pos- 
sible near the instability threshold and, for truly low- 
dimensional reduction, only for strongly confined sys- 
tems, with a discrete set of modes. The ST problem 
however is both unbounded and will operate in general 
far from threshold. On the other hand, since the growth 
is never saturated to a finite amplitude, any weakly non- 
linear analysis is necessarily limited to a rather early 
transient (57|. All the above techniques are thus of no 
much use for our purpose of studying the strongly non- 
linear dynamics of competing fingers in their way to the 
ST stationary solution. 

The basic point that we will exploit here to gain some 
analytical insight into the dynamics of the ST problems 
as a dynamical system is the fact that all exact solutions 
known explicitly for the idealized problem (do = 0) are 
defined in terms of ODE's for a finite number of parame- 
ters, and thus define finite-dimensional DS's in the phase 
space defined by those parameters |pq|. The complete 
ST problem, for any finite do, defines a DS in an infi- 
nite dimensional phase space. We will refer to this DS as 
S°°(do)- The limit do — ► defines a limiting DS which we 
will refer to as S°°(0 + ), which, as we will see, is different 
from S°°(0). 

The conformal mapping f(tu,t) of the reference unit 
disk in the complex cj-plane into the physical region oc- 
cupied by the viscous fluid z — x + iy = f(ui) has the 
form 



f(to, t) = — Iiilu + h(to, t) 



(13) 



where h(u>, t) is an analytic function in the whole unit 
disk, and therefore has a Taylor expansion 



h{u),t) = -Z%L ak(t)u> k 



(14) 



which is convergent in the whole unit disk. Inserting 
Eq.(|l4|) into the equation for the mapping f(co, t) we find 
an infinite set of equations of the form 



Ofc = gk(ao, —,ak]do). 



(15) 



In the co-moving frame of reference the precise form of 
the infinite set of equations fllq) is 
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Co = ^fov(0,t)d9 



and 
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(19) 



\vd u f(u>,t)\< 
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The infinite set of equations (15) defines the DS 
S°°(do). In the special case of strictly zero surface ten- 
sion, the DS S°°(0) can be explicitly solved for some 
classes of initial conditions. These classes define invari- 
ant manifolds of S°°(0) of finite dimension. In this con- 
text, finding explicit solutions implies identifying a spe- 
cific analytic structure of h(u>), with a finite number of 
parameters, which is preserved under the time evolution. 
If this condition is fulfilled, then a set of ODE's for those 
parameters can be closed, and defines a certain DS on 
a finite-dimension space. For instance, for do = the 
truncation of h(uj) into a polynomial form is preserved 
by the time evolution, so Eqs.(|l5|) themselves remain a 
finite set of ODE's. This simple case, however, is known 
to lead to finite-time singularities. The evolution is in 
general not defined after some finite time and cannot be 
considered as sufficiently well behaved as a DS's. On the 
other hand, classes of solutions have been reported which 
are smooth (non-singular) for all the time evolution. The 
corresponding conformal mapping takes the general form 

mm 
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h(u) = d(t) + Y^ Ij m (! - otj{t)w) 



(21) 



3=1 



where jj are constants of motion with the restriction 

J2 j=i lj — 2(1 — A) where A is the asymptotic filling- 
fraction of the channel occupied by fingers. If all 7j are 
real the evolution is free of finite-time singularities, and 
if any jj has an imaginary part then finite-time singu- 
larities may appear for some set of initial conditions (see 
Sect. VI C). Although this form of the mapping contains 



where 



all orders of the Taylor expansion of h(uS), it defines finite 
dimension invariant manifolds, since the superposition of 
logarithmic terms of Eq.(pij) is preserved under the dy- 
namics, this means that a closed set of ODE's for the 
finite number of parameters a,j(t) can be found. In ad- 
dition, the region which is physically meaningful is that 
in which \a.j\ < 1 (including the equal sign allows for 
the limiting case of infinite fingers, and makes the phase 
space compact). The DS defined by Eq.(pTl) in the 2N- 
dimensional hyper- volume will be denoted as L 2N ({jj}) 
For the sake of discussion throughout this paper it is 
important to have in mind that modifying parameters 
{7j}, which are constants of motion under the dynam- 
ics defined through Eq.(O) corresponds to varying ini- 
tial conditions in the phase space of S°°(0), while, from 



the standpoint of the finite-dimensional DS's denoted by 
L 2N ({jj}) it corresponds to changing the DS itself, that 
is, changing the ODE's obeyed by the dynamical vari- 
ables. In this sense, {7^} label a set of DS's defined on a 
2N- hyper- volume jay | < 1. 

Following Ref . [B8| , the key idea is to look for the sim- 
plest of the DS defined above which contains the three 
physically relevant fixed points, namely, the planar in- 
terface (PI), the single finger ST solution (1ST) and the 
double Saffman- Taylor finger solution (2ST) . We will call 
this minimal DS as L 2 (A, 0) or simply L 2 (X), since it 
has only one constant of motion, namely A. In Ref. |3 
we proposed to compare the global flow properties in 
phase space of such DS with those of a corresponding 
two-dimensional dynamical system defined by the reg- 
ularized problem. The latter was obtained by restrict- 
ing S°°(do) to a one-dimensional set of initial conditions 
properly chosen in such a way that the invariant manifold 
of 5°°(0) which defines L 2 (X) was tangent to S°°(0 + ) at 
the PI fixed point |i0|. The resulting DS S 2 (do) was then 
shown to have a topological structure nonequivalent to 
that of L 2 {\). In particular in the limiting case, S 2 (Q + ) 
intersected L 2 (l/2) not only at PI but also at the other 
fixed points 1ST and 2ST, and furthermore, at the two 
full trajectories connecting PI with 1ST and 2ST respec- 
tively. The basic conclusion was then that the regularized 
and the idealized problem were intrinsically different. 

In this paper we will see that the conclusions drawn 
from that comparison do hold beyond that 'minimal 
model' analysis. The extent and interpretation to which 
those conclusions must be understood will be more pre- 
cisely stated. We will also discuss on the implication of 
those results on the issue of dynamical selection. 



reducing the number of variables of the problem, but it 
is important to remark that it plays no role in prevent- 
ing finger competition, as seen in the regularized problem 
with identical symmetry. The case a'(t) = is equiva- 
lent to the time-dependent ST finger solution KM, and 
a"{t) — corresponds to the double time-dependent ST 
finger. For|a(i)| <C 1 the interface consists of a sinusoidal 
perturbation of the planar interface. d(t) affects the over- 
all interface position but does not affect its shape, and it 
turns out to be irrelevant for the present discussion. 
It is convenient to parameterize the phase space using 
] the variables u = 1-a" 2 and r = (a /2 +a" 2 -l)/(a" 2 -l). 
Then, the new phase space is the square [0,1] x [0,1]. 
With these variables (u, r) the Saffman- Taylor single fin- 
ger corresponds to the point (0, 1), the double Saffman- 
Taylor finger corresponds to (1,0) and the planar inter- 
face to (1,1). Substituting the ansatz Eq. ( f22[) in Eq. 
(O), we obtain that the temporal evolution in the new 
variables is given by 

/ , 3r — 4 — qr(l — ru) .„„. 

"^a-") l +9 T g (u,r) (23) 

, = 2r(1 _ r) 3r-2(l + ra ) +g (l-ru)(2-r) 



where 



l + gT g (u,r) 



1 



T g (u, r) = (1 - g)(2r + g(2r - 1)) - -(1 - gfru 

-gur 2 (l + g(ru~3)) (25) 

with g = 1 — 2 A = const. The detailed study of the two- 
dimensional dynamical system defined by these equations 
will be the object of the following section. 



IV. THE TWO-FINGER MINIMAL MODEL 



B. Study of the dynamical system 



A. The model 

The simplest exact time-dependent solution of Eq. 
(O) containing the three physically relevant fixed points: 
the planar interface (PI), the single Saffman- Taylor 
(1ST) fixed point and the double Saffman-Taylor (2ST) 
fixed point was introduced in Ref. Bq] and reads 

f(ui, t) = - lnw + d(t) + (1 - A) ln(l - a(t)u) 

+ (1- A)ln(l+a(i)*w) (22) 

where A is a real- valued constant in the interval [0,1], 
a(t) = a'(t) +ia"(t) and d(t) is real. The relevant phase 
space for a given A is the first quadrant of the unit circle 
in the (a', a") space. The other three quadrants describe 
interface configurations that are equal or symmetrical 
to the interfaces contained in the first quadrant. The 
interface described by this mapping consists generically 
of two unequal fingers, axisymmetric and without over- 
hangs. The axisymmctry of fingers simplifies the analysis 



The phase portrait of the dynamical system defined by 
Eqs. ( pffi - 125| ) with A = 1/2 is shown in Fig. 1. Its most 
interesting feature is the fact that the basin of attrac- 
tion of the Saffman-Taylor single finger is not the whole 
phase space. The separatrix between the basin of attrac- 
tion of the ST finger and the rest of the flow starts in 
the planar interface fixed point and ends in a new fixed 
point located at u* — and r* = 2A/(1 + A). The flow 
in the region below the separatrix is not attracted to a 
single fixed point, but it evolves to a continuum of fixed 
points, the line r = 0. All these fixed points correspond 
to stationary solutions of the zero surface tension prob- 
lem consisting of two unequal fingers advancing with the 
same velocity. The fingers have different width and tip 
position, the widths Aj.,2 given by 



Ai. 



1 ± — cotg 

7T 



1 



(26) 



and the tip separation in the propagation direction x be- 
ing 



A, = (1 - A) In 



l + yr 

1-vT 



(27) 



The complexity of the phase portrait shown in Fig. 1 is 
not at all evident with the original variables (a', a"), as 
can be seen in Fig. 7a. Had we not used the parame- 
terization (u, r), the non-trivial structure of the ST lim- 
iting solution would have remained hidden. The points 
along the line u = (equivalent to the point a' = and 
a" = 1 in the original variables) can be assimilated to 
the Saffman- Taylor finger because the dominant finger 
has the asymptotic ST finger shape for the correspond- 
ing A. All these points differ only in the evolution of the 
second finger, which always evolves to a needle (a finger 
of zero width), but of different length. When a trajectory 
approaches the 1ST fixed point from the two finger region 
the second finger does not disappear but remains frozen, 
that is, with lim^oo Ls/ ' Ll = where L$ and Ll are the 
lengths of the short and the long finger respectively. At 
the other extreme of the u = line, the point u = r = 0, 
the two fingers satisfy lim^oo Ls/ Ll = 1, but the dis- 
tance between the tips goes to infinity. The fixed point 
(u* , r*) corresponds to a new type of asymptotic station- 
ary solution of the zero surface tension problem, and it 
consists of two fingers with unequal positive velocities. 
Their length ratio satisfies lim^oo Ls/ Ll = 1/3 for any 
A. 



PI 




2ST 



FIG. 1. Phase portrait of the minimal model with A = §. 
The one-finger (resp. two-finger) region is above (below) the 
short-dashed line. For the region above the long-dashed line 
the secondary finger has zero growth rate while for the region 
below the secondary finger has finite growth rate. 

One important feature of the phase portrait, regarding 
the competition process, is the position of the line sepa- 



rating the regions where the small finger has zero (upper) 
or non-zero (lower) flux (see Fig. 1). A crossing from the 
region with non-zero flux to the one with zero flux would 
be the characteristic signal of 'successful' competition in 
the sense defined in Section II. But for the physically rel- 
evant value of A = 1/2, no trajectory crosses this line 
from the lower region or, equivalcntly, fingers with finite 
growth rate keep it finite for all time. Therefore this min- 
imal model does not exhibit successful competition, for 
A = 1/2. 

Although A = 1/2 is the physically selected value of 
the problem regularized with do — ► 0, the study of val- 
ues of A =/= 1/2 may be relevant to other purely Lapla- 
cian growth problems (without surface tension), such as 
DLA or needle-like growth. For fluid fingering problems 
it may also be relevant to situations where the selected 
finger width differs from 1/2 due to external perturba- 
tions (bubbles at the tip) or anisotropy. It is thus inter- 
esting to extend the analysis to arbitrary A. 



PI 




2ST 



FIG. 2. Phase portrait of the minimal model with A = yjj. 
The one-finger (resp. two-finger) region is above (below) the 
short-dashed line. For the region above the long-dashed line 
the secondary finger has zero growth rate while for the region 
below the secondary finger has finite growth rate. Note that 
there are trajectories crossing from below the line separating 
the zero and finite growth rate regions. 



The position of the fixed point r* that separates the 
flux depends monotonically on A, spanning the whole 
segment (0,1). For A > 1/2, r* approaches 1ST. Con- 
sequently the basin of attraction of 1ST is reduced. For 
A < 1/2 the behavior is the opposite: r* departs further 
from 1ST and the basin of attraction of the single fin- 
ger grows. In addition, as A decreases, a critical value 



A c = 1/3 is reached for which r* crosses the line sep- 
arating the zero and non-zero growth rate of the small 
finger. This crossing occurs at the point u = 0, r = 1/2. 
Therefore, for A < 1/3, r* will be located below that 
line, and some trajectories will cross the line from be- 
low: these trajectories exhibit successful competition, by 
definition. In Fig. 2 it is shown the phase portrait for 
A = 1/10 and trajectories crossing the line from below are 
depicted. Despite this fact, such successful competition 
is rather anecdotic, since it appears in a region far away 
from the neighborhood of the double finger, where the 
two fingers have a similar flux. The trajectories crossing 
the line from below correspond to a second finger with a 
small growth rate, thus the amount of flux totally elimi- 
nated is relatively small and quantitatively not significant 
in the competition process compared to the regularized 
problem. In general, what we can say is that under zero 




FIG. 3. Individual growth rates Aip\(t) and Aip^it) of the 
two fingers for the minimal model with (a) A = | and (b) 
A = yjj. For the (a) case the trajectory is attracted to a 
configuration of two unequal fingers and for the (b) case the 
trajectory reaches the Saffman- Taylor finger with successful 
competition. The growth rate of the secondary finger is the 
curve below. 



surface tension dynamics, narrow fingers present 
'stronger' competition than wide fingers. In Fig. 3 we 
show the evolution of the flux for a trajectory ending 
in a point of the continuum (case with A = 1/2) and 



a trajectory exhibiting successful competition (case with 
A = 1/10). 



C. Comparison with the regularized dynamics 

We are interested in the comparison between the do = 
dynamics and the do ^ one. The dynamical sys- 
tem defined by the mapping Eq. (K2) is referred to as 
L 2 (A). From now on we will restrict the analysis to 
most relevant case of A = 1/2. We introduce the follow- 
ing construction in order to compare the dynamics with 
and without surface tension: consider a one dimensional 
set of initial conditions (t = 0) of the dynamical sys- 
tem (|25J - p5|) surrounding the planar interface (PI) fixed 
point u = 1, r = 1. An example for this set is sim- 
ply a quarter of circle of small radius centered at (1,1), 
given by (1 - Rcos(6),l - Rsm(9)) with < 6 < tt/2 
and R -C 0. For a fixed A, the mapping Eq. (|2^) ap- 
plied to this set of initial conditions defines a continuous 
uni-parametric family of interfaces, where 8 is the pa- 
rameter in this example. The evolution up to t — > oo 
and backwards to i -> -oo of the case with do finite 
from this initial uni-parametric set spans a compact two- 
dimensional phase space (referred to as S 2 (do)) embed- 
ded in the infinite dimensional space of the problem with 
finite surface tension, S co (d a ). It is known from all ex- 
perimental and numerical evidence that for finite surface 
tension the subspace S 2 (do) must contain (at least) three 
fixed points. These three fixed points are the (unstable) 
planar interface (PI), the (stable) ST single finger (1ST') 
and a saddle fixed point that corresponds to the degen- 
erate double ST finger (2ST'), where the primes denote 
the finite surface tension case. The Saffman- Taylor fixed 
point (1ST') is known to be the universal attractor for fi- 
nite surface tension, so all the trajectories starting in the 
planar interface (PI) end up at 1ST'. The saddle fixed 
point 2ST' has an attracting manifold of lower dimen- 
sion, defined by a" = 0, and will govern the dynamics 
of finger competition. We define the space S' 2 (0 + ) as the 
limit of S 2 (d ) for d -» 0. 

In the neighborhood of the planar interface, namely the 
linear regime, the surface tension acts as a regular per- 
turbation and the regularized problem for small do con- 
verges regularly to the zero surface tension case. Thus, 
the manifolds L 2 (l/2) and 5 2 (0 + ) must be tangent at 
the PI fixed point u = 1, r = 1. Moreover, from selection 
theory we know that lST'^lST and 2ST'^2ST in the 
limit do — ► 0, provided that A = 1/2 is set for the zero 
surface tension manifold. Therefore, L 2 (l/2) and S 2 (0 + ) 
must be tangent at 1ST and 2ST. 

We have shown that the fixed points of the regular- 
ized space S 2 (0 + ) are contained in the zero surface ten- 
sion one L 2 (l/2), but the L 2 (l/2) contains additional 
fixed points that are not part of the regularized problem 
S 2 (0 + ). As we have seen in the previous section, these 
additional fixed points are the continuum of two-finger 



steady solutions. Therefore, it is important to stress that 
the flow for the two-dimensional subspace S 2 (0 + ) defined 
above is not topologically equivalent to the flow L 2 (l/2) 
of the minimal model. More precisely, the do — dynam- 
ical system contains a new fixed point that separates the 
phase space in two basic regions: the set of points end- 
ing up at 1ST and the set not ending up at f ST. Despite 
the fact that we do not know explicitly the trajectories of 
the regularized problem in the phase flow S 2 (0 + ), we can 
conclude that the two cases are fundamentally different. 

We have chosen the minimal model with the require- 
ment of including the double degenerate ST finger fixed 
point because it is essential to the competition process. 
In fact, from a physical point of view it is clear that 
the 2ST fixed point must govern the crossover from the 
growth regime to the competition one and its saddle point 
structure is necessary to the physical phenomenon of fin- 
ger competition, in the sense that it has an attracting 
manifold (containing the PI fixed point) associated to the 
growth regime, and an unstable manifold associated to 
the competition regime. However, the zero surface ten- 
sion dynamics fails to reproduce the saddle-point struc- 
ture necessary to describe the crossover from growth to 
competition. This failure is not a particularity of the 
minimal model but a generic property of the zero surface 
tension problem, as we will see in the following sections. 

The two-dimensional dynamical system (E3LE4LE5|) de- 
fines a structurall y u nstable flow, in application of the 
Peixoto theorem |56|], due to the existence of the con- 
tinuum of fixed points at the line r = 0. This structural 
instability is clearly related to the unphysical behavior of 
the minimal model, and its effects will be felt by solutions 
formally close to (P3), as it will be shown in the following 
section. According to the theory of dynamical systems, 
the introduction of an arbitrarily small perturbation to 
the dynamical system equations ( P3| , p4p5| ) could act as 
an unfolding of the problem, suppressing the continuum 
of fixed points and replacing the non-hyperbolic 2-equal- 
finger fixed point by an isolated one with a saddle-point 
(hyperbolic) structure. In this way, the unfolding of the 
phase portrait defined by the minimal class (22) would 



space is the following: 
f{u,t) = -lnw + d(t) 



have the same topology as the physical problem. There- 
fore, the natural step is to introduce a perturbation to 
the mapping ( p2| ) that should unfold the phase portrait 
and still be solvable. The next section is devoted to such 
perturbed ansatz. 



f (1 - A + ie)ln(l - a{t)w) 
-(l-A-ie)ln(l + a(*)*u;) (28) 



where e is a real positive and is a constant of motion. If 
e is set to zero then the minimal model Eq. (E3) is re- 
trieved. Solutions of this type have been studied before 
for instance in Refs. p9| , |62| , f33| . This mapping describes 
generically two unequal axisymmetric fingers, with the 
symmetry axis located in fixed channel positions sepa- 
rated a distance w, half the channel width. The main 
morphological difference between the interface described 
by the minimal class and the interface obtained from 
Eq. (|2g) is that the interface of the modified model may 
present overhangs. This can be understood from a geo- 
metrical point of view in the following manner: well de- 
veloped fingers are separated from each other by fjords of 
the viscous phase, and the width and orientation of these 
fjords is determined by the constant term I — A + it. If 
the constant term is real, i.e. e = 0, the centerline of 
the fjords is parallel to the channel walls and the fin- 
gers do not present overhangs, but if the constant term 
has an imaginary part (e ^ 0) then the fjords form a 
finite angle with the walls [Q . As a result of the incli- 
nation of the fjords the fingers may present overhangs. 
An example of these solutions is shown in Fig. 4, with a 
series of snapshots of the corresponding time evolution. 
The class of solutions Eq.(28) contains also the single 
finger Saff man- Taylor solution [a 1 = 0) but, remarkably 
enough, the introduction of a finite e has removed the 
double Saffman- Taylor finger solution. 




FIG. 4. Time evolution of a configuration with A = 1/2 
and e = 0.1. 



EXTENSION WITHIN TWO DIMENSIONS: 
SEARCHING FOR AN UNFOLDING 

A. Modified minimal model 



A possible modification of the ansatz (122) which is solv- 
able and preserves the two-dimensionality of the phase 



As usual, the constant A is the asymptotic width of the 
advancing finger. The natural phase space in this case 
is the unit circle, \a\ < 1, but we will restrict the study 
to a' > because the a' < region can be obtained by 
a 7r rotation of the a' > region. Physically, this rota- 
tion or the replacement a — > — a corresponds to a shift 
of the interface by an amount 7r (half the channel width) 
in the y direction. Thus, the semi-circle a' < contains 



the same interfacial configurations and dynamics than 
the a' > region after a trivial transformation. For the 
minimal model the zeros ojq of d u f(u, t) laid outside the 
unit circle, but for the modified minimal model Eq. (|28| ) 
the situation is different. For |a| < 1 a zero of d^f(u),t) 
can be inside the unit circle. The position of the zeros in 
this case is 



w 



-i(Xa" - ea') ± y/{2X - l)|a| 2 - (Xa" - ea'f 



(2A-l)|a| 2 



(29) 



O) 



for A 7^ 1/2 where a' — Re a and a" — Im a. For the 
particular value A = 1/2 the position of the zero is 



1 



UJ 



2i(Xa" - ea') ' 



(30) 



It can be shown that for any A and t / a wo can 
be found such that |a>o| < 1 for some |a| < 1. For in- 
stance, with A = 1/2 the curve |wo(a)| = 1 is the line 
a" = — 1 + 2ec/, which clearly intersects the unit circle 
\a\ — 1, enclosing a region where |wo| < 1- As a con- 
sequence of the presence of a zero inside the unit circle 
the parameter space \a\ = 1 contains unphysical regions, 
where the mapping Eq. (J2^) does not describe physically 
acceptable situations, with self-intersection of the inter- 
face associated to the fact that the mapping is not a single 
valued function. One of these regions is defined by the 
existence of a zero loq of 9 w /(w,i) inside the unit circle. 
In this region of phase space the interface crosses itself 
at one point, describing a single loop (see an example 
in Fig. 5). Most remarkably a second unphysical region 
containing interfaces with two intersections cannot be so 




FIG. 6. Time evolution of a configuration with a double 
crossing of the interface, with A = = an d e = \ ■ The solid line 
corresponds to t — with a = 0.85 + i0.4 and the short-long 
dashed line to t = 3.0. (The curves are plotted with its mean 
x position shifted arbitrarily for better visualization). 



Substituting the mapping Eq. (|28|) in Eq. ( |12[) it can 
be checked that this ansatz is a solution and that e is a 
constant preserved by the dynamics. The system of dif- 
ferential equations resulting from the substitution takes 
the form 

d + 41m[a(l-7)]{dlm[a]+lm[7«]}+ \a\ 2 {2X- 1) x 

{d\a\ 2 + 2 Re[ 7 a*d]} = 1 + |a| 4 + 4Im[a] 2 

2{d Im[a] + Im[ 7 d]} + 2d Im[a(l - 7)] + 

2 Im[a(l - l)]{d\a\ 2 + Re[7a*d]} + 

2(2A- l)\a\ 2 {dlm[a] +Im[7d]} = 4(1 + |a| 2 )Im[a] 

2Ad|a| 2 + 2Re[7a*d] = 2\a\ 2 

where the time-dependence of a(t) and d(t) has been 
dropped for sake of clarity and 7 = 1 — A + it. Eqs.(pl|- 
p3| ) can be integrated explicitely and the corresponding 
solutions for the variables d(t) and a(t) — a'(t) + ia"(t) 
take the form 



(31) 



(32) 
(33) 



(3 = d{t) - In a{t) + (1 - A - ie) ln(l - 
+(l-A + *e)ln(l 



t + C = Xd(t) + (1 - A)ln|a(i)| -earctan 



a(t)\ 2 ) 
a(t) 2 ) 
a"(t) 



a'(t) 



(34) 
(35) 



where C is a real-valued constant and (3 is a complex- 
valued constant. Notice that there is no apparent indi- 
cation of the pathological situations described above in 
the form the explicit solutions above. The study of the 
dynamical systems defined by Eqs. (f54|,[35|) is the object 
of the next section. 



FIG. 5. Interface with one crossing, with one zero inside 
the unit circle. 



B. Study of the dynamical system 



easily detected since, in this case, the zeros of 9^/(o;,i) 
lay outside the unit circle. Zero surface tension solutions 
displaying this feature were also reported in Ref. |63[ . 
Fig. 6 shows a configuration with this double crossing. 



The addition of an imaginary part ie to the constant 
(1 — A) modifies dramatically the phase portrait of the 
minimal model, as can be seen in Fig. 7, where the phase 
portraits for e = (now in the variables a', a") and 
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e = 0.1 are depicted. The phase portrait of the modified 
minimal model is qualitatively different from the e = 
one, as a direct consequence of the structural instability 
of the latter case [|36|, which implies that an arbitrary 
perturbation of the equations yields a flow which is not 
topologically equivalent (notice that a perturbation of an 
initial condition of the infinite-dimensional space of inter- 
face configurations is represented here as a perturbation 
of the equations themselves, that is, a displacement in the 
space of dynamical systems.) One could have expected 
that the introduction of this e would have provided an un- 
folding of the phase portrait of the minimal model into 
a structurally stable one (within the integrable class of 
mappings), hopefully with the saddle-point connection 
between the unstable and the stable fixed points, as cor- 
responds to the the physical case with do ^ 0. The phase 
portrait of the regularized flow (which is obviously the 
natural unfolding) would be similar to that of e = in 
Fig. 7a, except that all trajectories other than the line 
a" = would end up symmetrically to the upper ST 
fixed point or the lower one. Notice that in this rep- 
resentation, the 1ST fixed point has been split in two 
1ST(R) and 1ST(L), corresponding to whether the right 
or the left finger approaches the single finger attractor. 




1ST(L) -1 



FIG. 7. Phase portrait of the minimal model and the mod- 
ified minimal model. A = | for both plots, the regions to the 
right of the dotted lines correspond to two-finger configura- 
tions (a) e = 0; note the continuum of fixed points (marked 
with a thick line) on \a\ = 1. (b) e — 0.1; the straight line in 
the lower left corner is a line of finite-time singularities and 
the two fingers have equal length on the dashed line. 



These two solutions correspond to having the ST fin- 
ger located at two different positions (the symmetry axes 
of the fingers), owing to the translational invariance as- 
sociated to the periodic boundary conditions. Since a 
7r-shift in the transversal y-direction must yield an equiv- 
alent configuration, the identification of any point in the 
semicircle with the diametrically opposed which has re- 



duced the actual phase space in half, implies also that 
the 1ST(R) and 1ST(L) must be topologically identified 
as the same point. Approaching one or the other thus 
means approaching from the left or from the right. With 
this identification the contact with the problem with rigid 
boundaries is clearer, since there, the attractor is clearly 
unique but the flow must also be symmetrically split into 
two parts, corresponding to whether the left or the right 
finger wins, owing to the symmetry of the system under 
parity (see a more detailed discussion in section VI. D). 
Unfortunately, we must conclude that the resulting phase 
portrait for the modified minimal model does not provide 
the correct unfolding. This is particularly remarkable if 
one takes into account that, in two-dimensional systems, 
structurally stable dynamical systems are dense |56|| . On 
the contrary, the perturbed equations contain finite-time 
singularities and, although they remove the continuum 
of double-finger fixed points, they also miss the equal- 
finger fixed point, which is an essential ingredient of the 
regularized flow. 

In Fig. 8 we plot the phase portrait for e = 0.5 and 
the different regions of phase space. For any other e the 
flow is topologically equivalent but the shape and size of 
the different regions varies smoothly. The line of finite- 
time singularities collapses towards the lower fixed point 
1ST(L) in the limit e — > as shown in Fig. 7b. In the 
absence of the 2ST fixed point, the splitting of flow is 
made possible by the existence of the line of finite-time 



1ST(R) 1 




1ST(L) -1 



FIG. 8. (a) Phase portrait of the modified minimal model 
with A = 1/2 and e = 1/2. (b) Plot of different regions 
of phase space of case (a). The grey regions correspond to 
single finger interfaces and the other regions to two finger in- 
terfaces. Regions Ha and lib differ in which of the two fingers 
is larger. Regions III and IV are unphysical regions described 
in the text. The straight boundary of region III is a line of 
cusp singularities. 



singularities. Instead of a separatrix between the respec- 
tive basins of attraction of 1ST(R) and 1ST(L), there is 
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an intermediate, non-zero measure region, connected to 
the PI fixed point, whose evolution ends up at that sin- 
gularity line, defined by the condition |a>o| = 1- Similarly 
to the finite-time singularities occurring for polynomial 
mappings, this line is reached in a finite time and is as- 
sociated to the formation of a cusp at the interface. The 
evolution is not defined after that time. The flow in the 
region below the singularity line (region III of Fig. 8b), 
defined by \ojq\ < I, is actually well defined although 
it describes evolution of unphysical interfaces which in- 
tersect themselves forming a loop, (see Fig. 5). Their 
evolution originates and ends at different points of the 
singularity line. The region IV has double crossings of 
the interface (see Fig. 6) and also originates at the sin- 
gularity line but, remarkably enough, it evolves asymp- 
totically towards the ST finger despite their unphysical 
double crossing at the tail of the finger. This double- 
crossing is removed in a finite time in some subregion of 
IV and it remains up to infinite time in another subre- 
gion. Incidentally, this clearly illustrates how dangerous 
it may be to infer a physically correct dynamics from the 
fact that the interface evolves asymptotically towards a 
single ST finger, and that zero surface tension solutions 
must be dealt with extreme care since smooth and appar- 
ently physical interfaces may contain elements that yield 
them physically unacceptable when the time evolution is 
considered cither forward or backward. 




FIG. 9. Time evolution of a configuration with a double 
crossing of the interface, showing the dynamical removal of 
the crossing. In this case, A = i and e = ~. The initial in- 
terface with a = 0.865 + i0.2 is the solid line, and the dotted 
line is the interface at t = 0.5. The mean x position of the 
interface is shifted for better visualization. Time reversal of 
this evolution, corresponding to stably stratified Hele-Shaw 
flow, defines a finite-time interface pinch-off. 



The double-crossing removal in some of the above so- 
lutions has some implications in the general study of 
topological singularities associated to interface pinch-off 
in fluid systems (for a recent review see Ref. pm and, 
in the context of Hele-Shaw flows, see for instance Ref. 
p5[). Consider the stable Saffman- Taylor problem, in 
which the viscous fluid displaces the inviscid one. The 



planar interface is stable in this case and is the attrac- 
tor of the dynamics. The conformal mapping obeys then 
an equation formally equivalent to Eq.([12), that applies 
to the unstable Hele-Shaw flow, with the only difference 
that time is reversed, t is substituted by —t, in Eq.(fl2|). 
Therefore, the dynamics of the stable case is obtained 
from the unstable one simply by a time reversal. As a 
consequence, the double-crossing removal we observed in 
the original problem encompasses a prediction of a finite- 
time interface pinch-off in the stable configuration of the 
problem, for some class of initial conditions. A similar 
pinch-off phenomenon for zero surface tension dynamics 
was detected numerically by Baker, Siegel and Tanveer 
]63[ for other types of mapping singularities. Our result 
provides a very simple example of exactly solvable finite- 
time pinch-off. Notice that there is no singularity of the 
interface shape or velocity at the interface contact, so 
one could presume that surface tension may not affect 
significantly the phenomenon in this case, although this 
is an open question yet. 

The evolution of a trajectory ending up in the cusp 
line cannot be continued beyond the impact time to of 
the zero ujq with the unit circle \co\ = 1 because the cusp 
line attracts the flow from both sides, but the flow is in- 
deed well defined in the interior of region III, where we 
have -jet < in opposition to the physical region where 

—^ > 0, so, in a sense, the temporal evolution inside the 
singularity region is reversed. As a matter of fact, the 
graph a" (a') can be obtained and is smooth in all re- 
gions of phase space. Defining a — re l9 , its substitution 
in Eq. (fMj) yields after some algebra 



(10 
dr 



ArcosO (1 - A)(l- r 2 )sin# + e(l + r 2 )cos< 



. - r 2 1 + (2A - l)r 4 + 2Ar 2 cos 29 + 2er 2 sin 29 

(36) 

and from this expression the trajectory can be obtained 
also in region III. The fact that the modified minimal 
model does not yield an unfolding of the minimal one is 
more deeply stressed by the fact that the field of direc- 
tions defined by the above graph, even removing the sin- 
gularities through a proper time reparametrization and 
time reversal in region III, is still a structurally unstable 
flow. 

It is well known that the zero-surface tension problem 
is extremely sensitive to initial conditions: given a zero- 
surface tension solution at t = 0, another one can be 
found which is as close as desired to the interface of the 
first solution, and the evolution of the two solutions be 
completely different in general. This is a consequence of 
the ill-posedness of the initial- value problem |35[ , which 
is most clearly manifest in polynomial mappings, which 
may arbitrarily approximate any initial condition, possi- 
bly one which does not develop singularities, but them- 
selves always develop cusps. However, it is illustrative to 
see several striking examples of sensitivity to initial con- 
ditions within the class of logarithmic mappings which 
are much less predictive a priori . 
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Example 1 . Consider two initial conditions (a' 1: a'[) 
and (a'2,01'2) close to the PI fixed point, with |ai|, \a^\ <C 
1, which differ only in nonlinear orders of their mode 
amplitudes |66|. 




no signature of the different fate of the system could ap- 
parently be seen in the initial conditions. In both cases 
the extremely small differences associated to higher or- 
ders in the mode amplitudes have thus been crucial. The 
sensitivity to initial conditions of these examples is more 
striking for decreasing values of e, since the time in which 
the two evolutions stay close to each other increases as 
0(— lne). For instance, given an initial condition ao close 
to PI, the difference between the e = interface and the 
e — > one will remain of 0(e) for a time of 0(— lne). 
Later on in the evolution the differences between the two 
interfaces will be of O(l): the asymptotic shape of the 
e = case will be two unequal fingers while the shape of 
the e — *■ will be a single Saffman- Taylor finger. Simi- 
larly, for two initial conditions symmetrical to linear or- 
der such as in Example 2, with e — > 0, the differences 
between their interfaces will remain symmetric to 0(e) 
for a time of 0{— In e), but later they will lose symmetry 
and finally both will end up at the same fixed point, say 
the right one, even though one of the two evolutions has 
been favoring the other one, say the left one, for a long 
time (up to well developed fingers). 



FIG. 10. Evolution of two interfaces initially equal to 
linear order (see text), with A = 1/3 and e = 0.1. 
a(0) = 0.04619398 - iO. 01913417 for the solid line and 
a(0) = -0.04619398 -i0.00527598 for the dashed line. Upper 
left plot t — 0, upper right plot t — 2.0, lower left plot t — 4.0 
and lower right plot t = 6.0. 

One can easily choose (a^a'^) (with a^a'2 < 0, that is, 
considering not only the semi-circle a' > but the whole 
unit circle) such that the time evolution will be com- 
pletely different from the evolution of the original initial 
condition, even though the two initial conditions where 
equivalent to linear order. In Fig. 10 we show an explicit 
example. While the two initial conditions for the inter- 
face configuration cannot be distinguished in the scale of 
the plot, the final outcome is dramatically different. One 
of the evolutions is an example of successful competition, 
where the finger in the initial condition is eventually ap- 
proaching the ST solution, with a small secondary finger 
(not present in the initial condition) which is generated 
but screened out by the leading one. The other evolution 
is quite surprising since the secondary finger grows to the 
point of taking over and winning the competition. 

Example 2 . A similar situation is found if one com- 
pares two initial conditions equivalent to linear order up 
to a parity transformation. Pairs of initial conditions 
of this type, with the same values of A and e, can eas- 
ily be found within the same semicircular phase space, 
and since the dynamics is indeed symmetric under mir- 
ror reflection, one should not expect, in principle, a very 
different behavior, even though such points are not close 
to each other in phase space. Fig. 11 shows an example 
in which one of the evolutions is smooth, with a leading 
finger and a small one being generated, and the other 
generates a cusp in finite time. As in the first example, 




0.1)1) 



0,05 0,10 




FIG. 11. Evolution of two interfaces symmetric to linear or- 
der (see text), with A = \, e = 0.1, a(0) = 0.02724 -H0.03104 
for the solid line and a(0) = 0.02724 -i0.04193 for the dashed 
line. The unpper plot corresponds to i = and the lower to 
i = 4.19, when a cusp develops. 



Example 3 . In Fig. 12 we illustrate the effect of chang- 
ing e in initial conditions which are equivalent to linear 
order. Notice that in one case a cusp is generated at 
the secondary finger. In others the small fingers rapidly 
overcomes the large one, while for the smallest e the ini- 
tial finger seems to lead the competition. Remarkably, in 
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this last case the smaller finger will also take over after 
a much longer time. 




FIG. 12. Evolution of four interfaces equal to linear or- 
der, with A = f . e = 0.01 and q(0) = 0.04671 - i0.01291 
for the dotted line, the solid line is e — 0.1 and 
a(0) = 0.04619 - iO.01913, the short-dashed line is e = 0.3 
and a(0) = 0.04260 - J0.03137 and the long-dashed is e = 0.6 
and a(0) = 0.03472 - iO.04345. (a) t = 0, (b) t = 2.4, (c) 
t - 3.5, (d) t = 5.0. Note that at t - 2.4 the e = 0.6 interface 
develops a cusp. 

All the above examples have been chosen to empha- 
size the caution that is required when trying to use exact 
solutions to approximate the dynamics of the problem. 
A direct comparison of these solutions with numerical 
integration for very small surface tension would be re- 
quired in order to make a more quantitative assessment 
of the issue. This will be presented elsewhere |55| . In any 
case, it must also be stated that the class of logarithmic 
solutions does provide also qualitatively correct evolu- 
tions, not only of single finger configurations as stated 
in pa ], but also with two-finger configurations showing 
successful competition. An example of this is plotted in 
Fig. 4. Starting from the planar interface, during the lin- 
ear regime a bump starts to grow, followed generically 
by a second bump as the evolution enters the non-linear 
regime. The two fingers keep on growing for some time, 
until later on in the evolution, one of the fingers is dy- 
namically eliminated of the competition process and the 
other finger approaches asymptotically the ST finger so- 
lution. This general scenario is illustrated in figure 13a, 
where the individual growth rates of the two fingers A^i 
and Aip2 are plotted versus time, for two different initial 
conditions. 

For other initial conditions as generic as the previous 
one, however, the following phenomenon is observed: the 
small finger (with initially zero flux) of a configuration 
with two significantly unequal fingers increases its flux 
while the flux of the large finger decreases, until the flux 





FIG. 13. Individual growth rates Aipi(t) and Atp2(t) of the 
two fingers for the modified minimal model with A = § and 
e — 0.1, for two different initial conditions showing successful 
competition. For the (a) case the finger which initially has 
larger growth rate (and larger length too) wins the compe- 
tition. For the (b) case the finger which initially has lower 
growth rate (and lower length too) wins the competition, 
in opposition to the evolution with the regularized dynam- 
ics (small surface tension). 



of the initially small finger is higher than the flux of 
the other finger and finally the flux of the initially large 
finger reaches zero: it has been suppressed from the com- 
petition. This is opposite to what it would expected a 
priori from simple physical considerations. In fact, one 
would expect that for well developed fingers the larger 
one wins the competition, at least if the distance between 
the two finger tips is large before the process begins. This 
anomalous competition dynamics is illustrated in figure 
13b, where it can be seen that initially only one finger has 
a finite A^i, and it grows until a second finger develops 
and begins to grow, as indicated by the appearance of a 
non-zero Aip2 ■ The second finger grows together with its 
growth rate and surpasses the first one, which is eventu- 
ally suppressed from the competition process. This is in- 
dicated by Aipi going to zero. This example is important 
as a case where there is successful competition (finger co- 
alescence) to the Saffman- Taylor asymptotic solution but 
with a completely wrong dynamics. In fact it can be seen 
that the zero surface tension evolution departs from the 
regularized trajectory much before the small finger takes 
over the competition (through the impact of a daughter 
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singularity |59[). The winning finger with the regularized 
dynamics is thus the losing finger with the zero surface 
tension one. 

Again, in the limit e — > these phenomena appear 
even more dramatically, as a consequence of the struc- 
tural instability of the minimal model. In this limit, for 
a 0{— lne) time we will observe two unequal fully de- 
veloped fingers advancing with a fixed tip distance, but 
eventually the presence of finite e will 'activate' the com- 
petition process and one of the two fingers will reduce its 
growth rate until fully suppressed from the competition. 
If a" (0) > the suppressed finger will be the small one, 
but if a"(0) < the dynamically suppressed finger will 
be the large one. 



C. Comparison with the regularized dynamics. 

In order to compare the do = dynamics with the 
physical case of do ^ 0, as we have done for th e e = case, 
we use the construction introduced in Sect. 



IV C 



Con- 
sider a one dimensional set of initial conditions (t = 0) 
of the form Eq. (128) surrounding the planar interface 
(PI) fixed point a' — 0, a" — 0, for a fixed A and e. 
We choose the points of this set infinitesimally close to 
a = 0, and therefore the interface is in the linear regime. 
The do t^ time evolution from t = -oo to t = oo of 
this set spans a compact two dimensional phase space 
and defines a dynamical system S 2 (do), embedded in the 
infinite dimensional S°°(do). We define the dynamical 
system S' 2 (0 + ) as the limit of S 2 (do) for do — > 0, inter- 
secting the one dimensional set at t = 0. 

Since in the linear regime the regularized problem for 
vanishingly small do converges regularly to the do — so- 
lution, then the dynamical system L 2 (X, e) of do = and 
5 2 (0 + ) must be tangent at the PI fixed point, a = 0. Ac- 
cording to solvability theory the selected width A of the 
ST finger in the limit d -» is A = 1/2. Then, L 2 (l/2, e) 
and S 2 {0+) must intersect at 1ST(R), a' = 0, a" = 1 
and 1ST(L), a' = 0, a" = -1. In addition, from Rcf. 
p7{ it follows that the single-finger solution with d — > 
converges uniformly to the the time-dependent Saffman- 
Taylor finger solution with A = 1/2. This implies that 
L 2 (l/2,e) and 5 2 (0 + ) intersect not only in the fixed 
points that they have in common, but also on the line 
a' = (that corresponds to the time-dependent ST solu- 
tion). For the dynamical system S 2 (0 + ), the basins of at- 
traction of 1ST(R) and 1ST(L) are two-dimensional and 
finite, and therefore there must be at least one separatrix 
trajectory between the two basins. This separatrix must 
end at a saddle fixed point (which does not exist in the 
phase portrait of the do = solution). It is reasonable to 
assume that this fixed point is the double ST finger fixed 
point (2ST). Thus, the topology of the flow defined by 
the dynamical system with d = 0, L 2 (l/2, e) is different 
from the flow of the dynamical system do — > 0, S 2 (0 + ): 
the flow for the regularized problem contains a trajec- 



tory and a fixed point that it is not contained in the flow 
defined by the modified minimal model, the trajectory 
starting at the planar interface PI fixed point and ending 
up at the 2ST fixed point. The phase flow of the modified 
minimal model with do — is qualitatively different from 
the phase flow of the regularized problem, do — * 0, and 
therefore the solution Eq. ( pq ) is unphysical in a global 
sense, what is to say, when a sufficiently large set of ini- 
tial conditions (spanning evolutions towards 1ST(R) and 
1ST(L)) is considered simultaneously. 

It is useful to define a projection of the trajectories 
from the regularized two-dimensional space 5* 2 (0 + ) onto 
the do = L 2 (l/2,e) space in order to make possible 
not only the comparison between the topology of the two 
flows but also the comparison between individual trajec- 
tories. It is obvious that there are many different pos- 
sible definitions, and the projection of the regularized 
flow 5 2 (0 + ) onto L 2 (l/2,e) is in that sense not unique, 
but the following discussion should not depend on the 
particular definition of the projection. One possible defi- 
nition of a useful projection is the following [p5[ : given a 
point £ of S 2 (0 + ) and a point a of L 2 (l/2, e), we define a 
'distance' D{ (a) between £ and a, that is, a distance be- 
tween an interface within 5* 2 (0 + ) and an interface within 
L 2 (l/2,e). One possible choice of Dq(cx) is the area en- 
closed or trapped between the two interfaces. Then, we 
define the projection a^ of C (that corresponds to an in- 
terface solution of the problem with do 7^ 0) onto the 
do = phase space L 2 (l/2, e) as the value of a that min- 
imizes the 'distance' D^(a) between a and £, or equiv- 
alently, that minimizes the 'distance' between the two 
interfaces, with the restriction that the average position 
of the two interfaces is the same to ensure that the pro- 
jection satisfies mass conservation. Then, once we have 
defined a projection a^ of a point of 5 2 (0 + ) it is straight- 
forward to obtain the projection of a particular trajectory 
£(£) in S 2 (0 + ), applying the definition of a^ to any point 
of the trajectory. This projection of a trajectory will be 
denoted by a^ (t) . 

The projection on the do = phase portrait of the 
do — ► separatrix trajectory as e -> will approach 
the time-dependent 2ST trajectory for e = 0, the line 



a" = 0, but a broad set of initial conditions of e 







located below the projection of the separatrix (that is, 
with a" < 0) will be attracted by the 1ST(R) fixed point 
under the do = dynamics instead of being attracted by 
the 1ST(L) fixed point as their position with respect the 
separatrix implies |55| . Thus, the evolution under the 
do = dynamics of this set of initial conditions will be 
not only quantitatively but also qualitatively incorrect, 
in the sense that the finger that wins the competition is 
the one that in the regularized dynamics loses. An ex- 
ample of this was given in Fig. 13b. More details will be 
presented elsewhere |p9[ . 

The motivation for the introduction of a finite ie term 
to the minimal model Eq. (22) was to obtain an un- 



folding of its non-hyperbolic fixed point structure, but in 
this section we have shown that the introduction of ie has 
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failed to unfold the phase portrait to the expected physi- 
cal topology, the saddle-point structure of the regularized 
problem. It has dramatically changed the topology of the 
flow obtained for e = 0, but the flow for e 7^ does not 
have the expected structuraly stable flow of the physical 
problem: an unstable fixed point, two stable fixed points 
and one saddle fixed point. Instead of this, the evolution 
of Eq. (Pq) with e / presents finite-time singulari- 
ties for a non-zero measure set of initial conditions. This 
can be understood as a consequence of the absence of the 
2ST saddle point, which controls the competition regime. 
Without this fixed point the separatrix trajectory neces- 
sary to separate the basins of attraction of ST(L) and 
ST(R) is not present and the only possible way to split 
the flow is through the existence of finite time singulari- 
ties. This is not a particularity of the mapping Eq. (Eq) 
but a more general feature of do =0 solutions. Below 
we will prove that, within the N-logarithms class, finite e 
implies finite-time singularities in the evolution of a non- 
zero measure set of initial conditions (see Sect. VIC). 
Besides the existence of finite-time singularities we have 
seen that, unlike the case e = solutions exhibiting suc- 
cessful competition are possible with e 7^ for A = 1/2. 
However, part of those evolutions are unphysical in the 
sense the winning finger may differ from the one with the 
regularized dynamics. 



VI. GENERALIZATION TO HIGHER 
DIMENSIONS 

This section is devoted to the study of solutions that 
define a dynamical system of dimensionality greater than 
two. We will show that the discussion of previous sections 
is not peculiar of dimension 2 but is extendable to higher 
dimensions. 



A. Non-axisymmetric fingers 

The solutions that have been studied in the previous 
sections, Eq. ( |22] ) and Eq. (|2q), have two pole- like singu- 
larities Wi,2 located at u>i = l/a and 102 — —1/a*. The 
property uj\ — — cj^ reduces the dimensionality of the dy- 
namical system to two and also forces the axisymmctry 
of the fingers. If the singularities 0^2 are not related, 
then the phase space has one additional dimension and 
the fingers are not axisymmetric. 

We will now study the ansatz 



f(uj,t) = -I11L0 + d(t) + (1 - A) ln(l - c*i (t )u>) 

(l-A)ln(l-a a (t)w) 



(37) 



where otj{t) = ct'Jt) + ia"(t) are two complex quantities 
satisfying |aj| < 1. This ansatz is an exact solution of 
Eq. ( p"2| ) and it can be proved that the evolution is free 
of finite time singularities. 



From the evolution equations for 0^2 in the case 
A = 1/2 it is found that the dynamics satisfy 



= Const. 



(38) 



£ -1 U.n 



£ -1 \.jLc\ 



Writing a in polar coordinates, aj 
straint reduces to the simpler form 

0i + 9 2 = C. 



r,-e J , this con- 



(39) 



The constant C depends on the initial condition, but its 
value can be fixed arbitrarily using the property of rota- 
tional invariance of the ansatz Eq. (|37|) : the transforma- 
tion a.\ — > aie 1 ^ and a.2 — > o^e 11 ^ is equivalent to a trans- 
lation of the interface a distance </> in the y axis direction. 
For the minimal class studied previously (Eq. (E3)) the 
value of C was C = it. The existence of the constraint 
Eq. ( p9| ) reduces the dimensionality of the problem from 
four to three, with variables n, r% and 9\ — #2, implying 
that the dynamical system is actually three-dimensional. 
In order to study the dynamical system defined by the 
substitution of the ansatz Eq. (pTJ) into the evolution 
equation ( J12| ) we introduce the variables z = r exp i9 = 
ai + «2 and p = ol\ol2- Moreover, we use the arbitrari- 
ness of C to choose the simpler value C = 0. Then, the 
dynamical system in the variables (r, 9, p) reads 



4r 



1 - r 2 + p(2 + r 2 ) cos{29) - 3p 2 
4_ r 2 

2(l-p 2 )+r 2 (pcos(2(9)-l) 



p = Ap 

9 = -2psm(29). 



(40) 



With these variables the axisymmetric case studied in 
previous sections is contained as a the particular case 
9 = 0. From the last equation it is immediate to see that 
9 decreases monotonically with time (since, for C = 0, 
p > 0), and tends to zero for long times, 9^0, converg- 
ing to the axisymmetric case 9 = 0. Therefore, we can 
conclude that the axisymmetric case is attracting the dy- 
namics of the non-axisymmetric one and, asymptotically, 
its dynamics reduces to the axisymmetric one. 

Therefore, we have shown that the dynamics of the 
minimal model is not a particularity of a zero-measure 
subset but the general behavior of the non-axisymmetric 
neighborhood of the original axisymmetric class. This 
neighborhood has a non-zero measure within the three 
dimensional system. Similarly, it can be shown that the 
basin of attraction of the Saffman- Taylor finger for the 
ansatz Eq. ( |37| ) is also relatively small. 

It is worth noting that non-axisymmetric fingers can 
be obtained also in 2d. A conformal mapping describing 
non-axisymmetric fingers is 

f(u, t) = - \nui + d(t) + {l-\+p + ie) ln(l - a(t)u>) 

+(l-\-p-ie)ln(l + a*(t)u>) (41) 
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where < p < 1 — A. However, the lack of axisym- 
metry caused by the introduction of p does not change 
qualitatively the phase portraits obtained for p = e = 0, 
the continuum of fixed points present for e = is not 
removed by the introduction of a finite p and the finite- 
time singularities that appear for £ ^ are also present 
when p ^ 0. Therefore, the results obtained for 2d in the 
previous sections are not modified at all if we relax the 
condition of finger axisymmetry. 



B. Perturbations which change finger widths 



2ft we 



The second type of modification of the ansatz 
have studied is the following. Consider 



f(uj,t) = -\nuj + d(t) + (l- A)ln(l-ai(t)w) 

+(l-A)ln(l-a a (*)a>) (42) 
+2 (A- \ s )ki(l -5(t)u) 



with initial conditions ai(0) 



;(0) - a(0), < 



A,A S < 1 and |<5(0)| < 1. Note that for 5(0) = this 
ansatz reduces to the minimal model (p2[). From substi- 
tution of this ansatz into the evolution equation (112) it is 
obtained that Eq. (|42| ) is a solution with A and A s con- 
stants. From the dynamical equations it can be proved 
that the asymptotic configuration of this ansatz consists 
of one or two fingers, with asymptotic filling fraction 
equal to A s . But if \S(0)\ <C |ck(0)| then the interface will 
be initially almost identical to the one obtained within 
the class (p2|) with the same a(0) and A, and its evolution 
will remain close to the one obtained with the minimal 
class for a time that will increase with decreasing |<5(0)|. 
Therefore, given a small enough |<5(0)|, starting from the 
planar interface a configuration with one or two fingers 
(depending on the initial conditions) of total width A will 
develop. Later on, as \5\ grows and approaches 1, the to- 
tal width will change from A to X s for long enough time. 
The ansatz ( f42| ) thus describes an interface that changes 
the filling fraction of the fingers from A to A s . The same 
phenomenon will appear with any other of the solutions 
described in this paper (and in general in pole-like solu- 
tions) if a term of the type 2(A— X s ) ln(l— 5(t)u) is added, 
and in particular it will appear in the single finger config- 
urations obtained setting a(0) and 5(0) purely imaginary. 
Note that in this case if 5(0) = the ansatz describes 
exactly the time-dependent Saffman- Taylor finger. This 
changing-width phenomenon of do = solutions has been 
known for long [|48| , but it has been recently claimed |2q | 
to be responsible of the known width selection observed 
with finite surface tension both experimentally and nu- 
merically. The idea was that, although solutions of arbi- 
trary A exist in the absence of surface tension, these are 
unstable under some perturbations which trigger the evo- 
lution towards the A = 1/2 solution. Since the present 
paper is basically emphasizing the unphysical dynamics 
of the idealized (do = 0) problem, in direct contradiction 



with Ref. pq] , we feel compelled to briefly comment on 
this respect here. The basic argument of Ref. J2£| is as 
follows, in terms of the parameterization of the interface 
used by the author: a term of the form i/j,<p in the confor- 
mal mapping is always unstable under the substitution 
ificf) — ► fi\n(e l< ^ — e). The introduction of such perturba- 
tion then leads to the a = case, which corresponds to 
A = 1/2. In Refs. |29| , |30|| it was pointed out that, with 
the same degree of generality, equivalent perturbations 
exist which lead to any desired A, and therefore the con- 
clusion that A = 1/2 is the only attractor is incorrect. 
It is argued |31|] that the latter class of perturbations is 
different form the former since they increase the num- 
ber of logarithmic terms in the conformal mapping and 
therefore modify the dimension of the subspace of solu- 
tions. This objection is somewhat misleading since such 
partitioning of classes of solutions in terms of the number 
of logarithms is arbitrary and not intrinsic. This can be 
seen by choosing a different reference region to confor- 
mally map the physical fluid. Instead of mapping it into 
the semi- infinite strip |2q ], the mapping into the interior 
of the unit circle avoids the confusion on the dimension 
of the subspace of solutions. Thus, the perturbation pro- 
posed in Ref. [28J is equivalent to choosing A s = 1/2 in 
the ansatz ([42 ), but it is manifest in this formulation 
that there is nothing special with this particular choice 
of A s . Perturbations leading to any finger width A s occur 
with the same genericity. Therefore the instability of the 
point 5 = is not related to the steady state selection 
phenomenon. 



C. Finite-time singularities within N- logarithms 
solutions 

In this section we will prove that any solution of the 
N-logarithm class [OT that does not have only real con- 
stant parameters presents finite time singularities, that 
is, it contains a non-zero measure set of initial conditions 
which develop singularities at finite time. 

Consider a conformal mapping function f(to,t) 

f(u>, t) = - In u + d(t) + (Ai + ie) ln(l - a\ (t)u) (43) 
+(A 2 - ie) ln(l -a 2 (t)uj) 

where Ai + A2 = 2(1 — A), e > and ai,2 are com- 
plex with |cki 2 1 < 1- The mapping f(ui,t) must satisfy 
d u f(w, t) ^ for \uj\ < 1. If any zero ivq of d w f(w, t) hits 
the unit circle \u)\ = 1 then the interface develops a cusp. 



For the ansatz (M) d u f(u>,t) reads 

1 (Ai + ie)ai _ (A 2 - ie)a 2 
lu 1 — a.\U) 1 — a 2 0J 

Thus, the position of the zero luq of 9oj/(woji) is 

_ — (Ai+ie — 1)qi— (A2— «e — l)a2 
^0 - 2Q!Q 2 (2A-1) 



(44) 



± 



■ v /((Ai+^-l)ai+(A2-^-l)a2) 2 -4aia 2 (2A~ , . 

2aia 2 (2A-l) V 4t V 
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If, for some value of 01,2, |«i,2| < 1, the zero uiq is in- 
side the unit circle, then the ansatz (|4^) will present fi- 
nite time singularities for some sets of initial conditions. 
Therefore, if \ujq\ < 1 the interface will develop a cusp. 
Setting ai,2 = ae iSl - 2 and 2 -6 1 = -28 with 5 < 1 the 
position of the zero (keeping up to linear terms in 8) is: 



U) 



, A±(1-A) 
q(2A-1) 



iSe 



-i0 2 



a(2A-l) 



A 2 - 1 - ie ± 



A - 1 + A(A 2 - ie) 



1-A 



0(8 2 ) 



and the modulus of the minus solution (the one with 
smaller modulus) reads 



N> 



=p 



1-A 



0(<5 2 )]. 



(47) 



In consequence, for a close to 1 we obtain |u>o| < 1, one 
of the zeros is inside the unit circle in a finite neighbor- 
hood of Qi = «2 = e l6 ■ Thus, the mapping ( |43| ) presents 
finite time singularities for some initial conditions inde- 
pendently of the value of e and A12, and the measure of 
this set is non-zero. 

Now we consider a generic mapping with N > 2 loga- 
rithmic terms of the form: 



N 



/(«,*) 



lnw 



d(t) + $>iMl-«i(*)w) 



(48) 



j=i 



where 7^ = Aj + ilj are constants of motion with the 

restriction X);=i7i = 2(1 — A). If we choose ctj = a\ 
for 1 < j < k and a,- == 0:2 for fc + 1 < j < N we 
recover the mapping ([43|). Therefore, the N-logarithm 
solution ( (48[ ) contains initial conditions that develop a 
cusp with this subset of a>j, but the dimension of this 
subset is lower than the dimension of the phase space, 
implying that the measure of the set is zero compared to 
the whole phase space. To prove that the subset that de- 
velops a cusp is finite we choose now the following values 
for uf aj = exx 4- T]j for 1 < j < k and ctj = a 2 + r\j for 
k + 1 < j < N, with \r)j\ < 1, where \uj \ < 1 if r]j = 0. 
The equation d u f(ui,t) = reads 



1 k 



7i(«i + »?i) 



(ai + Vj)v 



N 



7j(«2+%) 

- (a 2 +%) w 



= 0. (49) 



This equation (EJ) reduces to Eq. (Esj) if all r/j = 
and it has N zeros if ?7j 7^ 0. Defining g(u> ) = <9 w /(u;,£) 
for ?7j = and G(u,ff) — d u f(u,t) for 77^ 7^ then 
G(u>,fj) = #(w) + 8G(uj,ff) where |5G(a>,?y)| < AT^ for 
|w| < i?, with AT and i? constants, and <?(wo) = 0. One 
zero u' of G(w, 77) can be written lj' = ujq + 8uj, and as- 
suming \8lo\ < C\ff\ with C constant the substitution of 
lo' in G(u>, fj) = yields 



3(^o) 



9w 



<W + £G(wo,?7) = 0. 



(50) 



The position of the zero is then: 

8G(uj Q ,ff) 



U) Q = U) 



() 



~duj 



where -£■ 



(51) 



is m- 



7^ 0. Therefore, the zero Wq of Eq. 

side a ball of radius o(\ff\) centered in u>q. If |wo| < 1, then 
choosing \ff\ small enough the zero will satisfy \oj' \ < 1: in 

neighborhood of (01,02) at least one zero of d u f(oj, t) 
is inside the unit circle, and the dimension of this neigh- 
borhood will be the same of the phase space. So we can 
conclude that any mapping of the form ( fig ) presents fi- 
nite time singularities for some sets of initial conditions 
of non-zero measure, provided that at least one pair of 
7j has a non zero imaginary part. 

Thus, the requirement that a mapping function of the 
form (Eq) is free of finite time singularities for any ini- 



tial condition ctj(0) is fulfilled if and only if Imfyj 



0, 



j = 1,...,N. But this restriction implies [B2| that for a 
wide range of initial conditions the asymptotic configura- 
tion is a N-finger interface with unequal fingers advancing 
at a constant speed, a situation fully analogous to the one 
discussed in Sect. IV. Then, if a mapping of the form 
( pq ) with Im[7j] = is chosen, then the dynamical sys- 
tem L 2N (jj) will have non- hyperbolic fixed-points (con- 
tinua of fixed points) and will lack the saddle-point struc- 
ture of the regularized problem. In order to completely 
remove the continua of fixed points it is necessary to set 
Im[7j] 7^ f62fl , but in this case we will encounter finite- 
time singularities and the saddle-point structure will not 
be present anyway. 

To sum up, we have shown that the features of the min- 
imal model and its extensions that make them unphysical 
(in a global sense) are not specific of their low dimension- 
ality. The features that make the solutions studied in pre- 
vious sections ineligible as a physical description of low 
surface tension dynamics for a sufficiently large class of 
initial conditions, are also present within the much more 
general N-logarithm family of solutions, and the conclu- 
sions drawn in previous sections can be generalized to 
that class. 



D. Rigid- wall boundary conditions 



It is worth stressing here that the use of periodic 
boundary conditions throughout this study, as opposed 
to the physically more natural rigid- wall boundary condi- 
tions, is not essential to the basic discussion. In connec- 
tion with the discussion of multifmger steady solutions, 
this point was raised in Ref. [p7[ and addressed in Ref. 
fl68[ . Here we will just recall that the choice of periodic 
boundary conditions is not only the simplest in terms of 
symmetry and dimensionality, but it is the relevant one if 
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one is interested in general mechanisms of finger compe- 
tition in finger arrays. In this sense, the study of the two- 
finger configurations in this paper refers to an alternating 
mode of two-finger periodicity in an infinite array of fin- 
gers, in the spirit of Ref. [p9|. For finite size-systems one 
can also argue that rigid-wall boundary conditions are 
included as a particular case of periodic boundary con- 
ditions in an enlarged system. That is, a channel with 
width W with rigid walls in mathematically equivalent 
to a channel of width 2W with periodic boundary con- 
ditions where auxiliary channel of with W is constructed 
as the mirror image of the physical one. The competition 
of two fingers in a channel with rigid walls at a distance 
W is in practice equivalent to a four-finger problem with 
periodic boundary conditions in a double channel. 

The only subtle point which we would like to stress 
is the apparent degeneracy of the single-finger attractor 
into a left ST finger and a right ST finger, as already 
pointed out in Sect. VB, and the possible relevance of 
this fact in connection with the saddle-point structure of 
the phase space flow. This degeneracy is inherited from 
the trivial continuous degeneracy associated to transla- 
tion invariance in the transversal direction, when peri- 
odic boundary conditions are assumed. In fact an ar- 
bitrary shift in the transversal direction yields a physi- 
cally equivalent configuration. When an initial condition 
is fixed, such continuous degeneracy is broken into two 
discrete spatial positions which are separated a distance 
W/2. The whole dynamical system is then invariant un- 
der translations of W/2. This is the reason why we only 
plotted a half of the disk in the phase portraits of section 
M Technically, the resulting dynamical system must be 
defined 'modulo- W/2', that is, identifying any configura- 
tion with the resulting of a W/2 shift. In the phase space 
defined by the variables (a' , a") one should identify any 
point with the resulting of a 7r rotation. In this way the 
two single-finger attractors do correspond to the same 
fixed point. With this identification, the ST finger is not 
degenerate and the flow becomes topologically equiva- 
lent to the corresponding one in a channel with rigid-wall 
boundary conditions. The two-finger configurations have 
thus the same structure, regardless of the type of side- 
wall boundary conditions. The flow starts at the PI fixed 
point and ends up at the 1ST fixed point. Between them 
there is a saddle point corresponding to the 2ST fixed 
point. This separates the flow in two equivalent regions, 
namely 'from the left' and 'from the right' of the saddle 
point. With zero-surface tension, the case of rigid walls 
exhibits the same problems, namely the occurrence of 
a (nontrivial) continuum degeneracy of multifinger solu- 
tions, and the existence of finite-time singularities. The 
important point we want to stress is thus that all the gen- 
eral conclusions drawn in this paper are valid if rigid-wall 
boundary conditions are considered. 



VII. DYNAMICAL SOLVABILITY. GENERAL 
DISCUSSION 



A. The physics of zero-surface tension 

The role of the zero surface tension solutions in the 
description of the dynamics of the nonzero but vanish- 
ingly small surface tension problem is now clearer. The 
do = dynamics is in general incorrect in a global sense, 
even if we choose solutions with the asymptotic width A 
given by selection theory. However, they have an impor- 
tant place in the description of the physical problem. It 
has been proved in Refs. [B5h37| that the solutions with 
do = converge to the do —* during a time 0(1 ), before 
the impact with the unit circle of the so-called daughter 
singularity at time td- In practice this implies that the 
do = dynamics is not only correct in the linear regime 
(where do acts as regular perturbation) but also quite 
deep into the nonlinear regime. After td nothing can be 
said a priori: as we have shown in the present paper, 
there are regions of the d = phase space correspond- 
ing to smooth interfaces that are physically wrong, but 
other regions are a good description of the evolution with 
finite (but very small) surface tension. For instance, in 
the neighborhood of the time-dependent Saffman- Taylor 
finger (the line a' = in the solutions (E2), (pq)) the 
do = evolution is qualitatively correct for finite sur- 
face tension, and even quantitatively correct in the limit 
do ^ (for A = 1/2). However, a question remains 
open: given a do — evolution smooth for all time and 
consistent with the results of selection theory, is it the 
limit of a do — > evolution? This question can be ex- 
plored numerically and is the subject of a forthcoming 
paper [ p5[ . Generally speaking the conclusion is that ex- 
act solutions including evolution of two different fingers 
which are compatible with MS theory, that is, evolving 
to a single finger with the width predicted by selection 
theory, and which do not exhibit any kind of singularity 
in the interface shape, may be dramatically affected by 
surface tension. The outcome of the competition, that 
is, which one of the two competing fingers will survive at 
the end, when an infinitesimally small surface tension in 
introduced may be the opposite of that of the zero sur- 
face tension case. This may happen in situations where 
fingers are significantly different from each other and is 
not an instability of a particular trajectory, but a generic 
behavior in a finite (non-zero measure) range of initial 
conditions within the integrable class. For that region of 
phase space, it is clear that the dynamics of finger com- 
petition is completely wrong for the class of integrable so- 
lutions. Nevertheless, there is also a class of initial condi- 
tions which have a qualitatively correct evolution includ- 
ing 'successful' finger competition in the sense defined in 
sections above (this possibility was incorrectly excluded 
in Ref. |38|], where the analysis was based on e = 0). 
Although strict convergence of the regularized solution 
to the idealized one may not occur in these cases, the 
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quantitative differences may be moderately small. Ac- 
tual convergence of some type can only be expected at 
most when there is only one finger along the complete 
time evolution. In summary, according to this scenario 
there are basically four classes of initial conditions within 
the most general integrable solutions, once those a priori 
incompatible with selection theory are excluded, namely 
(i) finite-time singularities forward or backward (or both) 
in time; (ii) asymptotically correct ST finger with wrong 
dynamics (the incorrect finger wins); (iii) asymptotically 
correct ST finger with qualitatively correct evolution (the 
correct finger wins although shapes may differ during a 
transient); and (iv) (unphysical) evolution towards mul- 
tifinger fixed points. It has to be added that, all of the 
above solutions plus those which are incompatible with 
selection theory are qualitatively and quantitatively cor- 
rect in the limit of small surface tension, until a time of 
order one which is always in the deeply nonlinear regime. 

As a general consideration it is worth remarking that 
fingers emerging from the instability of the planar inter- 
face when this is subject to noise are necessarily in the 
range of dimensionless surface tension of order one. A 
simple way to argue this point is that it is precisely sur- 
face tension which selects the size of the emerging fingers, 
since the fastest growing mode is that in which both sta- 
bilizing and destabilizing forces are of the same order. In 
these cases, surface tension is felt necessarily in the lin- 
ear regime, and the usefulness of the zero surface tension 
solutions in the early stages of the evolution is obviously 
more limited. 

Finally, from a physical point of view it is appropri- 
ate to recall that the presence of noise does modify the 
general picture of the fingering dynamics in the limit of 
small surface tension, as pointed out in Ref. Q. Al- 
though the ST finger is the universal attractor of the 
problem, the linear basin of attraction decreases with di- 
mensionless surface tension. In practice this implies that 
the interface approaches the ST finger but when it gets 
too close, noise triggers its nonlinear instability and the 
interface makes a long excursion (typically a tip splitting) 
before approaching again the ST finger. The considera- 
tions made in this paper concerning the limit of small 
surface tension thus imply that noise must be taken as 
sufficiently small. A careful discussion of the effects of 
noise, particularly in numerical simulation of very small 
surface tension will be presented elsewhere [pal. 



B. A Dynamical Solvability Scenario 

In Ref. [f38| we pointed out for the first time the dy- 
namical implications of the MS analysis when extended 
to multifinger fixed points. The idea of the Dynamical 
Solvability Scenario (DSS) was already latent in that dis- 
cussion. We pursued this extension of the steady state 
selection problem explicitly in Ref. p|,f59| , where we found 
that, in direct analogy to the single-finger case, the in- 



troduction of surface tension did select a discrete set of 
multifinger stationary states, in general with coexisting 
unequal fingers. Here we would like to discuss in what 
sense that analysis does provide a Dynamic Solvability 
Scenario. 

Before doing that, let us briefly consider an alterna- 
tive view of a possible DSS proposed by Sarkissian and 
Levine J32|. In Ref. |32|] , it was explicitly discussed with 
examples that exact solutions of the zero-surface tension 
problem did behave differently from numerical integra- 
tion of the small surface tension problem. At the end, 
the authors speculated with the possibility that surface 
tension could play a selective role in the sense that it 
could basically pick up the physically correct evolutions 
out of the complete set of solutions without surface ten- 
sion, in direct analogy with the introduction of a small 
surface tension selecting a unique finger width out of the 
continuum of stationary solutions. Since the class of non- 
singular integrable solutions is indeed vast and infinite- 
dimensional, it is not unreasonable to expect that one 
could approximate any particular evolution with finite 
surface tension with one of those solutions for all time. 
However, as recently pointed out in Ref. [Q, there is 
no simple way to determine which of those solutions is 
selected by any macroscopic construction. Furthermore, 
even if this were possible, one should still face the rather 
uncomfortable fact that the base of solutions defined by 
the superposition of logarithmic terms in the mapping, 
would itself correspond to unphysical (nonselected) so- 
lutions, as we have seen throughout this paper. Indeed, 
an initial condition defined exactly by a finite number 
of logarithms would have to be replaced in general by a 
solution with an infinite number of logarithms as the 'se- 
lected' solution which the (small) finite surface tension 
system tracks. 

From a more general point of view, a dynamical selec- 
tion principle understood as 'selection of trajectories' has 
an important shortcoming when considered within the 
perspective of a broader class of interfacial pattern form- 
ing systems. In fact, the solvability theory of steady state 
selection has turned into a general principle because its 
applicability to a large variety of systems, most remark- 
ably in the context of dendritic solidification || BJO)- 
However, it is only for Laplacian growth problems that 
exact time-dependent solutions are known explicitly, so 
there would be no hope to extend the above DSS as a 
general principle to those other problems. 

The DSS we propose here has a weaker form but it 
is susceptible of generalization to other interfacial pat- 
tern forming systems. The basic idea can be best ex- 
pressed in similar words to those recently used by Gol- 
lub and Langer |fl|] to describe solvability theory in a 
general context. They have nicely synthesized the sin- 
gular role of surface tension in the language of dynam- 
ical systems as to 'whether or not there exists a stable 
fixed point' |ll]]. In this context, our DSS extends the 
(static) solvability scenario in the sense that the singu- 
lar role of surface tension is precisely to guarantee the 
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existence of multifmger fixed points with a saddle-point 
(hyperbolic) structure. We have seen that the contin- 
uum of multifmger fixed points is directly related to a 
nonhyperbolic structure of the equal-finger fixed points. 
They imply directions in phase space were the flow is 
marginal, and this is so to all derivative orders. While 
in the traditional solvability scenario the introduction of 
surface tension does isolate a stable fixed point (a con- 
tinuum of single-finger fixed points turns into a stable 
one and a discrete set of unstable ones), now it isolates 
multifmger saddle points out of continua of multifmger 
solutions, as discussed in Ref. |39],|| (a continuum of n- 
finger fixed points turns into a hyperbolic fixed point with 
stable and unstable directions, and a discrete set of un- 
stable ones) . Since the saddle fixed points are defined by 
the degenerate n-equal-finger solutions, the stable direc- 
tions of the saddle-point are directly related to the stable 
directions of the single- finger fixed point, while the un- 
stable directions correspond to all perturbations which 
break the n-periodicity of the equal-finger solution. The 
most important stable and unstable directions, however, 
are those depicted in the two-dimensional phase portraits 
discussed in the above section, namely the 'growth' direc- 
tion connecting the planar interface and the n-finger fixed 
point, and the 'competition' direction connecting the n- 
finger fixed point to the single-finger fixed point J70| . No- 
tice that arrays of fingers emerging from the morpholog- 
ical instability of the planar interface are relatively close 
of n-periodic solutions as long as the noise in the initial 
conditions is weak and white, which guarantees that the 
most unstable (fastest growing) mode dominates in the 
early nonlinear regime. In these conditions, the system 
feels the attraction to the corresponding n-equal finger 
fixed point. This stage is what we called 'growth'. When 
the fingers are relatively large they start to feel the devi- 
ations form exact periodicity and start the 'competition' 
process. 

Note that, despite the formal analogy to the single- 
finger solvability theory, the reference to a the restoring 
of multifmger hyperbolicity by surface tension as dynami- 
cal solvability scenario is fully justified. Indeed, the local 
structure of the multifmger fixed point has a dramatic 
impact on the global (topological) structure of the phase 
space flow, as we have seen in simple examples. The 
existence of a small but finite surface tension thus deter- 
mines a global flow structure and it is in this sense that 
it 'selects' the dynamics of the system. 

The possibility of extension of this analysis to other 
interfacial pattern forming problems relies on the exis- 
tence of a continuum of unequal multifmger stationary 
solutions with zero surface tension. The fact that in the 
ST case the existence of those can be associated to a 
simple relationship between screening due to relative tip 
position and relative finger width (that is, a slower areal 
growth rate of the screened finger is compensated by its 
smaller width, resulting in an equal tip velocity), one is 
tempted to conjecture that similar classes of solutions 
must exist in other problems, for instance in the growth 



of needle crystals in the channel geometry || . Although 
this point should be more carefully addressed, it seems 
reasonable to expect that a DSS as presented above could 
be generalizable, to some extent, to other physical sys- 
tems. 



VIII. SUMMARY AND CONCLUSIONS 



We have developed a Dynamical Systems approach to 
study the dynamics of the Saffman- Taylor problem, bas- 
ing the analysis on the zero surface tension solutions. A 
minimal model has been analyzed, and from its phase 
flow we have concluded that it is unphysical. A detailed 
study of a perturbation of the minimal model within two 
dimensions has yielded the same conclusion. The un- 
physical behavior of zero surface tension solutions is a 
consequence of the non-hyperbolicity of the multifmger 
fixed points of the finite-dimensional dynamical system 
that they define, opposed to the saddle point structure 
of the regularized problem. Perturbations of the minimal 
model to higher dimensions confirm the generality of the 
conclusions reached in two-dimensional models. We have 
proved that the N-logarithms class of solutions presents 
finite-time singularities if the continua of fixed points are 
totally absent. From the analysis of zero surface ten- 
sion solutions we conclude that they are unphysical in a 
global sense, when sufficiently large classes of initial con- 
ditions are considered simultaneously, because they lack 
the correct topology of the physical flow, structured in 
terms of a saddle-point connection between the unstable 
and the stable fixed points. This does not exclude that, 
for some sets of initial conditions, the zero surface ten- 
sion dynamics might be correct, not only qualitatively 
but even quantitatively, but it is not possible in practice 
to know it a priori by any simple means. We have illus- 
trated with several examples that although the asymp- 
totic behavior may be correct (evolution towards a single 
ST finger) the intermediate dynamics may be completely 
wrong, or even physically meaningless, such as for the 
existence of interface crossings. We have also illustrated 
the sensitivity to initial conditions when approximating 
physically relevant situations with different integrable so- 
lutions. As a by-product we have also obtained some 
practical results concerning zero-surface tension dynam- 
ics which may be relevant to Laplacian growth problems, 
for instance in relation to the interplay of screening ef- 
fects and finger widths. We have introduced precise def- 
initions of 'growth' and 'competition'. With the proper 
definition 'successful' competition, we can state for in- 
stance that, in the absence of surface tension, narrow 
fingers do compete more efficiently than wide ones. We 
have also found explicit solutions which lead to finite- 
time interface pinch-off in the stable configuration of the 
problem. 

The detailed comparison of the dynamics with zero 
and non-zero but very small surface tension requires a 
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careful numerical study and can be analyzed in terms of 
the daughter singularities formalism developed in Refs. 
|pq , p7| . As a matter of fact it can be shown that the 
zero surface tension problem and the vanishingly small 
surface tension regularization differ dramatically even in 
regions where the former is nonsingular, in the sense that 
non-zero measure regions of phase space have a different 
outcome of the competition (namely, which one of two 
competing fingers survives) in the two cases. A detailed 
study of this point will be presented elsewhere |5g] . 

Finally, we propose a Dynamical Solvability Scenario 
relevant in principle not only for viscous fingering prob- 
lems but also of applicability to other pattern forming 
problems. Within this DSS the role of surface tension 
as a singular perturbation is to isolate multifinger saddle 
points out of the continua of multifinger fixed points, as 
shown previously in Ref. [§9]||. This extends the tradi- 
tional solvability theory applied to steady state selection, 
where surface tension did also isolate a unique (stable) 
hyperbolic fixed point out of a continuum of nonhyper- 
bolic ones. In that case the isolated fixed point was the 
global attractor of the problem. In the present exten- 
sion, the introduction of surface tension does isolate a 
unique n-equal finger fixed point out of each continuum 
of n-finger fixed points, with both stable and unstable 
directions. By restoring this saddle point local structure 
the topology of the phase space flow is modified, so the 
introduction of surface tension has a deep impact on the 
global phase-space structure of the dynamics. It is in this 
sense that this scenario can be considered as a dynamical 
solvability theory. 
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